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Measurements of the fractional momentum loss (Sioss = Spt/pt) of high-transverse-momentum- 
identified hadrons in heavy ion collisions are presented. Using 7r° in Au-|-Au and Cu-fCu collisions 
~ 62.4 and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy 
Ion Collider and and charged hadrons in Pb-fPb collisions measured by the ALICE experiment at 
the Large Hadron Collider, we studied the scaling properties of S'loss as a function of a number of 
variables; the number of participants, A^part, the number of quark participants, Aqp, the charged- 
particle density, dNch/dp, and the Bjorken energy density times the equilibration time, eBjTo. We 
hnd that the pr, where Sioss has its maximum, varies both with centrality and collision energy. 
Above the maximum, Sioss tends to follow a power-law function with all four scaling variables. The 
data at = 200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common 

scaling of Sioss with dNoh/dp and EBjTo, lending insight on the physics of parton energy loss. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

It has been firmly established that in relativistic heavy 
ion collisions a hot, dense medium is rapidly formed, ca¬ 
pable of interacting with the high pr partons produced 
in primordial hard scattering and makingthem lose some 
energy while traversing the medium [iH^. Such ener^ 
loss in the medium was first predicted in early I980’s Q- 
Quantifying this energy loss is an important issue, be¬ 
cause it is directly connected to the properties of the 
medium. However, this is not straightforward since nei¬ 
ther the original parton energy, nor that of the decel¬ 
erated one is easily accessible. Back-to-back photon-jet 
pairs in principle give access to both the initial and fi¬ 
nal parton energy, but such events are rare, because they 
are suppressed by a factor a, the electromagnetic cou¬ 
pling constant. Measurement of jets give more complete 
information on the parton energy loss, however, their 
measurement is challenging, particularly at high multi¬ 
plicities and low parton pp. To circumvent this, high 
Pt hadrons are often used as proxies for jets (“leading 
hadrons”), and the parton energy loss in principle can be 
calculated by proper comparison of the invariant yields 
of hadrons in p+p and A+A at a given pp. For this pur¬ 
pose the p+p yields are usually scaled up by the expected 
number of binary nucleon-nucleon collisions in A+A, es¬ 
timated from a Glauber Monte-Carlo model, and in the 
absence of any initial or final state nuclear effects they are 
expected to coincide with the A+A yields. The partons 
have steeply falling momentum spectra, so if partons lose 
energy, that results in a shift of the momentum spectra, 
and the yield at a given pT will become suppressed Q. 
Utilizing this fact, the nuclear-modification factor (Raa) 
has become a widely used characterization of the energy 
loss which is defined as: 


Raa{pt) = 


{l/Nl^^)d^N^Jdppdy 

(Taa) X d^cT^Jp/dprdy ’ 


( 1 ) 


* Deceased 

^ PHENIX Co-Spokesperson: morrison@bnl.gov 
^ PHENIX Co-Spokesperson: jamie.nagle@colorado.edu 


where is the production cross section of the respective 
hadron in p+p collisions, {Taa) = (-^coii) /o’p'Jf* is the 
nuclear overlap function averaged over the relevant range 
of impact parameters, and (iVcou) is the number of binary 
nucleon-nucleon collisions computed with If Raa 

is unity, it is usually assumed that the yield measured 
in A+A collisions is explained by the primordial hard 
production as observed in p+p collisions with no nuclear 
or medium effect. If Raa < 1 (suppression) the A+A 
yield at a given pp is less than that expected from the 
scaled p+p. 

While the parton energy loss is expected to depend 
both on system size and collision energy, it is remarkable 
that Raa is very similar from = 62.4 to 200 GeV 

at the Relativistic Heavy Ion Gollider (RHIG) and up to 
2.76 TeV at the Large Hadron Gollider (LHG). The rea¬ 
son is that while the energy loss increases with increas¬ 
ing ^ysNN which would tend to decrease Raa , the power 
n in the Pt~'^ shaped spectra decreases (n = 10.6 for 
62.4 GeV l|, n = 8.06 for 200 GeV Au-kAu and n rv 6.0 
for 2.76 TeV 0) and provides a countervailing effect. A 
numerical calculation showed that the fractional energy 
loss of partons, AE/E, is indeed significantly different 
between LHG and RHIG even though the Raa is simi¬ 
lar Q. 

Instead of Raa one can employ the fractional momen¬ 
tum loss (S'loss) of high Pp hadrons as a measure of parton 
energy loss which should reflect the average fractional en¬ 
ergy loss of the initial partons {{AE/E) ^ S'loss)- Sioss is 
defined as 

pp _ AA 

Sioss = Spp/P p = —— pp ^ (2) 

Pp 

where p^^ is the pp of the A+A measurement and pi^ is 
that of the p+p measurement scaled by the nuclear over¬ 
lap function Taa of the corresponding A+A centrality 
class at the same yield of the A+A measurement. We 
calculate Sioss as a function of the original momentum of 
partons that are represented by p^. 

Under the assumptions that iVcoii scaling is applica¬ 
ble and fragmentation functions are unchanged from p+p 
collisions, Spp can be directly measured as the shift in pp 
needed to get the same yield {dN/dppdy) in A+A as the 
scaled p+p. 
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The PHENIX experiment published a study of the 
energy loss of partons by converting azimuthal angle 
((/))-dependent Raa with respect to the event plane to 
S'loss assuming that the spectra follow a power-law func¬ 
tion 0- That study found that S'loss scales with Lg, 
the distance from the center to the edge of the colli¬ 
sion area which the partons traverse, for all centrality 
classes for 3< pt <8 GeV/c, and also with the density- 
weighted path length pL /pcent where Pcent is the density 
at the center of the collision zone and the p is the den¬ 
sity at the given coordinate. The dependence of Sioss on 
centrality was also reasonably approximated by iVpart^^^- 
A similar study has been performed using Pb-|-Pb data 
available at LHC and Au-I-Au data from RHIC [ll|. The 
authors found that the scaling in [13 does not hold at pr 
higher than 10 GeV/c. Other recent publications tried to 
obtain ^-integrated 5'ioss without assuming the spectral 
shape [3,13 ■ It was found that S'loss varies by a factor of 
six from 62.4 GeV Au-|-Au to 2.76 TeV Pb-|-Pb collisions. 

These studies showed that the fractional momentum 
loss Sioss bas a major advantage over Raa, in that it 
allows for a direct comparison of parton energy loss be¬ 
tween different colliding systems and energies, because it 
eliminates the bias owing to the -variation of the 

exponent, n, in the power-law spectra of high pT parti¬ 
cles. 

These scaling studies are not a replacement for full 
quantum-chromodynamics calculations of parton energy 
loss that must include different quark and gluon admix¬ 
tures and their different fragmentation functions, initial 
state effects such as nuclear modified parton distribu¬ 
tion functions, and potentially modified harmonization 
effects. That said, since S'loss is merely a new represen¬ 
tation of the experimental measurements, any such the¬ 
oretical calculation would need to describe the observed 
scalings at the precision of the uncertainties. 

In this paper, we extend the previous studies of (p- 
integrated Sioss by including additional data sets both 
from RHIG and LHC and by plotting the fractional mo¬ 
mentum loss against several scaling variables to charac¬ 
terize the energy loss mechanism. We average over the 
event plane dependence to simplify the analysis. Sec¬ 
tion m describes the method of calculating Sioss and in¬ 
troduces the global scaling variables. In section UlIAl we 
present values for Sioss as a function of centrality for a 
variety of systems and energies. Section IHI HI presents 
the main result of this paper, which is the study of the 
scaling behavior of Sioss- We conclude in section ITVl 


II. DATASET AND ANALYSIS 

In this section we describe how fractional momentum 
loss is calculated and define the various scaling variables. 
A summary of the data is given in Table HI For RHIC 
energies, data from the PHENIX experiment for 7r° in 
Au-I-Au and Cu-|-Cu collisions both at = 200 GeV 

and 62.4 GeV were used [3. fl^ - fisjl ■ while for the LHC, 


TABLE I. Summary of data sets used in this analysis. The 
= 62.4 and 200 GeV data are from PHENIX at RHIC 
and the ~ 2.76 TeV data from from ALICE at the 

LHC. 


System 

particle 

y/SNN 

year 

pT range 

ref. 

Au+Au 

7r° 

200 GeV 

2004 

1.0-20 GeV/c 

[12] 

Au-I-Au 


200 GeV 

2007 

5.0-20 GeV/c 

[8j 

Cu-hCu 

7V° 

200 GeV 

2005 

1.0-18 GeV/c 

[13] 

P+P 


200 GeV 

2005 

0.5-20 GeV/c 

M 

Au+Au 

7V° 

62.4 GeV 

2010 

1.0-10 GeV/c 

m 

Cu+Cu 


62.4 GeV 

2005 

1.0-8.0 GeV/c 

[13] 

P+P 

TfO 

62.4 GeV 

2006 

0.5-7.0 GeV/c 


Pb+Pb 

h+^- 

2.76 TeV 

2010 

0.2-50 GeV/c 

[1^ 

Pb+Pb 


2.76 TeV 

2010-2011 

2.0-20 GeV/c 

[17] 

Pb+Pb 

7r° 

2.76 TeV 

2010 

0.5-11 GeV/c 

[18] 

P+P 

h+^- 

2.76 TeV 

2009-2011 

0.2-50 GeV/c 

[19] 

P+P 


2.76 TeV 

2010-2011 

2.0-20 GeV/c 

[17] 

P+P 

7r° 

2.76 TeV 

2011 

0.5-11 GeV/c 

[18] 


data on charged hadrons and pions in Pb-|-Pb collisions, 
both at = 2.76 TeV, measured by the ALICE ex¬ 

periment [la - H^ were used. To calculate the fractional 
momentum loss, p+p data are also needed: RHIC data 
were taken from [3 [13 , while LHC data were taken 

from [13. 


A. Fractional momentum loss 

Figured] shows the method of calculating the S'loss us¬ 
ing measured A+A and p+p spectra at the same collision 
energy. First, the , h^/~) cross section inp+p is 

scaled by Taa corresponding to the centrality selection of 
the A+A data. Second, the scaled p+p cross section is fit 
with a power-law function. Third, the scaled p+p point, 
pf^, corresponding to the yield at the Au-f Au point of in¬ 
terest, is found using the fit to interpolate between scaled 
p-hp points. The dpT is calculated as p^^-To obtain 
S'loss, the 6pT is divided by p^. 

It is important to realize that the effective fractional 
energy loss, Sioss, estimated from the shift in the pr spec¬ 
trum, is actually less than the real average energy loss at 
a given px- This is true because, for a given observed 
Px^, the events at much larger px with larger energy 
loss are lost under the events at smaller px with a corre¬ 
spondingly smaller energy loss owing to the steeply falling 
spectrum. We evaluated this bias to the S'loss measure¬ 
ment with a simple Monte Carlo calculation using the 
power of the spectra obtained in the measurements, and 
found that it is ~10% for collisions at = 200 GeV 

and 62.4 GeV, and -18% for = 2.76 TeV. This 

systematic effect is not reflected in the final data uncer- 
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taken from . The number of quark-participants is cal¬ 
culated for all systems as part of this work, as explained 
below. 

A Monte-Carlo-Glauber (MC-Glauber) model calcula¬ 
tion is used to obtain estimates for the number of 
nucleon participants at each centrality using the proce¬ 
dure described in [2^ . A similar procedure can be used to 
estimate the number of quark participants, A^qp, at each 
centrality [1^. The MC-Glauber calculation is modified 
such that the fundamental interactions are quark-quark 
rather than nucleon-nucleon collisions. The nuclei are 
assembled by distributing the centers of the nucleons ac¬ 
cording to a Woods-Saxon distribution. Once a nucleus 
is assembled, three quarks are then distributed around 
the center of each nucleon. In our model, we assume the 
spatial distribution of the quarks follows an exponential 
charge distribution as measured in electron-proton elastic 
scattering: 




(3) 


FIG. 1. (Color online) Method of calculating the fractional 
momentum loss (S'loss =5pt/pt)- This plot is for illustration 
only; uncertainties are not shown. The procedure: (1) scale 
the p+p data by Taa corresponding to the centrality selection 
of A+A data, (2) fit the p+p data and choose the scaled p+p 
point closest in yield to the A+A along the fit,(3) calculate 
the difference of scaled p+p and A+A transverse momenta, 
SpT = p^ — Pt^, at the same yield. 


where a = \/vllrm = 4.27 fm“^ and Vm = 0.81 fm is the 
rms charge radius of the proton [2^ . The coordinates of 
the two colliding nuclei are shifted at random relative to 
each other by a vector 6, the impact parameter, which 
covers an area larger than the maximum possible impact 
parameter. A pair of quarks, one from each nucleus, 
interact with each other if their distance d in the plane 
transverse to the beam axis satisfies the condition 


tainties. 

The uncertainties of the S'loss are obtained as follows. 
We first estimated the errors of yields for the A+A and 
the p+p points in three categories; the quadratic sum of 
the statistical and pr-independent systematic uncertain¬ 
ties (“Type A”), py-correlated systematic uncertainties 
(“Type B”), and the overall scale uncertainties which al¬ 
low all the data points to move to the same direction with 
a certain fraction of the central values (“Type C”). The 
Type B is the quadratic sum of the systematic uncertain¬ 
ties related to the measurement of 7r° for the PHENIX 
result, including those of photon identification efficiency, 
energy scale, and background subtraction. The Type G 
is the quadratic sum of the Taa and p+p normalization 
uncertainties in this analysis. The uncertainties for the 
A+A and p+p points in three categories are separately 
summed in quadrature, and projected to the p^ axis us¬ 
ing the p+p fit function. 


B. Number of Nucleon and Quark Participants 

To study the systematics of fractional momentum loss, 
we introduce several scaling variables. Here we briefly 
describe how the number of nucleon participants (A^part) 
and quark participants (A^qp) [1^ are obtained. The 
A^part for the Pb-|-Pb collisions at = 2.76 TeV was 


d< 



(4) 


where tr™®* is the inelastic quark-quark cross section, 
which is varied for the case of nucleon-nucleon collisions 
until the known inelastic nucleon-nucleon cross section is 
reproduced; this cr™®' is then used for the A-l-A calcula¬ 
tions. The inelastic quark-quark cross sections are tabu¬ 
lated in Table [TTl Figure [5^ shows the number of quark 
participants as a function of the number of nucleon par¬ 
ticipants [ 2 ^. The relationship is nonlinear, especially 
for low values of A^part- The nonlinearity is clearly seen 
in Fig.[5}D where the ratio of the number of quark partici¬ 
pants to the number of nucleon participants as a function 
of the number of nucleon participants is shown. 


TABLE II. The inelastic quark-quark cross sections used 
for each collision energy to reproduce the inelastic nucleon- 
nucleon cross section. 


(GeV) 

(mb) 

(mb) 

2760 

64.0 

18.4 

200 

42.3 

9.36 

62.4 

36.0 

7.08 
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FIG. 2. (Color online) (a) The number of quark participants as a function of the nnmber of nucleon participants. The error bars 
represent the systematic uncertainty estimate on the MC-Glauber calculation. The dashed line is a linear fit to the 200 GeV 
Au+Au points with Apart > 100 to illustrate the nonlinearity of the correlation at low values of Apart, (b) The ratio of the 
number of quark participants to the number of nucleon participants as a function of the number of nucleon participants. The 
error bands represent the systematic uncertainty estimate on the MC-Glauber calculation. This figure is reproduced from [2^. 


C. Charged Particle Multiplicity 


D. Bjorken Energy Density 


Another scaling variable used is charged particle mul¬ 
tiplicity, or multiplicity density, dNch/drj, measured at 
midrapidity (y ~ ~ 0). This quantity is closely related 

to the gluon density, dNgiuon/dy [1^, as well as to the 
number of participating nucleons Apart, which in turn is a 
measure of the system size. In a previous publication 
it has been shown that 


dAch/dy OC Apart“ (5) 


where a=1.16 in Au-|-Au collisions at = 200 GeV. 

For the RHIC data dNrh/dri values were taken from the 
PHENIX experiment [^, [^, where charged particle 
multiplicities are measured in the \r]\ < 0.35 pseudora¬ 
pidity region in two pad chamber detectors [26j in zero 
magnetic field. For the LHC data dNrh /d? 7 . values are 
quoted from the ALICE publication [^, where charged 
particles are measured in their silicon-pixel detector and 
quoted in the restricted |ry| < 0.5 pseudorapidity range. 


Finally, we introduce a measure of the energy density. 
In relativistic heavy ion collisions, the Bjorken energy 
density is frequently used for this purpose [l^. The 
Bjorken energy density is defined as 


1 dEj' 
roA_L dy 


( 6 ) 


where tq is the proper time when the QGP is equilibrated, 
A± is the transverse area of the system. The A± can be 
written as ^ cr^ay, where and <7y are the widths of 
X and y position distributions of the participating nu¬ 
cleons in the transverse plane, and was estimated using 
a Monte-Carlo Glauber simulation The equilibra¬ 
tion time To is strongly mo del-dependent, therefore, we 
decided to use Ebjto as a scaling variable, which then con¬ 
tains only well-established experimental quantities. The 
measured dET/dy is converted to dEx/dy by applying 
a factor that compensates the phase space difference be¬ 
tween rapidity and pseudorapidity which is obtained by 
a simple numerical calculation. The factor is found to be 
1.25 for = 62.4 GeV and = 200 GeV [H, 

and 1.09 for = 2.76 TeV The uncertainties 

on these scale numbers are ^3%. The dEx/dr} for the 









= 2.76 TeV Pb+Pb collisions are obtained from 
the literature [^ . 


III. RESULTS AND DISCUSSION 


The numerical values of the scaling variables defined 
in the previous section are listed in Table Ell 

a. 0-4 
i-" 0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 
0 


4 6 8 10 12 14 16 18 20 

pj" (GeV/c) 

FIG. 3. (Color online) dependence of Sioaa for 7r° in 
200 GeV Au+Au collisions from (solid symbols) 2007 data @] 
and (open symbols) 2004 data from the PHENIX experiment 
at RHIC for pr <10 GeV/c [T2|. The error boxes corre¬ 
sponding to Type-B errors are not shown for Year-2004 data, 
but the magnitude are same as the ones for Year-2007 data. 
Saya{TAA © pp uorm) are Type-C errors and show the absolute 
amount that the data points would move. 


+ PHENIX 
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A. pt dependence of the fractional momentum loss 

Figure [3] shows the px dependence of the fractional 
momentum loss of 7r° for various centralities in Au+Au 
200 GeV collisions, using 2007 data Q. The error bars 
represent the projection of Type A uncertainties to the 
axis, while the boxes are the same projection of Type 
B uncertainties. SgysiTAA © PP norm) shown in the fol¬ 
lowing plots stands for the projection of Type C uncer¬ 
tainties to the Px axis. Note that Ssys{TAA © PP norm) 
indicate the absolute amount that the data points would 
move. 

The 2007 data set has been analyzed only above px 
= 5 GeV/c, which also limits the px where S\oss can be 
extracted. For lower px the 2004 data were used 
and the results are shown in open symbols in Fig. [3] The 
consistency of Raa from 2004 and 2007 data has already 
been shown in Fig. 11 of [H- The same consistency can 
be seen in the extracted iSioss- In the central collisions 
S'loss is slightly increasing up to GeV/c, then flattens 


out and finally decreases at the highest measured px- As 
expected, S'loss increases monotonically with centrality. 

We show the fractional momentum loss of 7r° for vari¬ 
ous centralities in Gu+Gu 200 GeV collisions in Fig. S) 



pP" (GeV/c) 

FIG. 4. (Color online) p^ dependence of Sioss for tt® in 
200 GeV Cu+Cu collisions using the spectra measured by 
PHENIX at RHIC in 2005 [l^. 3sya{TAA © pp norm) are 
Type-C errors and show the absolute amount that the data 
points would move. 

We already found in a previous publication that Raa 
is similar at the same Apart between Gu+Cu and Au+Au 
collisions at = 200 GeV [l^. The Apart for 0%- 

10% centrality in Cu+Gu collisions is similar to the one 
for 30%-40% centrality in Au+Au collisions. We can see 
that the S'loss is similar in these collision from Figs. [3] and 

SI 

The fraction of hard-scattering is smaller and therefore 
results in a steeper px spectrum at = 62.4 GeV. 

Figure [5] shows the fractional momentum loss of 7r° for 
various centralities in Au+Au 62.4 GeV collisions. 

The Sioss is much smaller than at 200 GeV even for 
the most central collisions. Note that soft production in 
A + A collisions still contributes to the p^ range of 2- 
6 GeV/c, where Raa is not reaching to its minimum Q . 
In the Sioss, this will result in smaller values. Figure [B] 
shows the Sioss of 7r° for various centralities in 62.4 GeV 
Cu+Gu collisions Q- 

The trends are similar for the Cu+Cu and Au+Au 
collision data. Note that in the 62.4 GeV data set the 
systematic uncertainties from 7r° reconstruction, overall 
energy scale and trigger efficiency were larger [l3l| than 
in the 200 GeV Au+Au data, which explains the larger 
overall systematic uncertainties. It is again interesting 
to mention that within the uncertainties, the 0%-10% 
Cu+Cu collisions give the similar Sioss as the 20%-40% 
Au+Au collisions even at this energy. 

In Fig. El we show the fractional momentum loss 
for charged hadrons in Pb+Pb collisions at ,/s „ = 

2.76 TeV measured by the ALICE experiment (l^. fl^. 
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TABLE III. Global variables for Au+Au and Cu+Cu collisions at RHIC from PHENIX 0,1,El [3 and Pb+Pb collisions at 
the LHC from ALICE Ei,E3,|33]- 


Collision 

\/^NN 

Centrality 

Apart 

Aqp 

dN^^/dr) 

esj-ro [GeV/fm^] 

Au-I-Au 

200 GeV 

0%-5% 

353±10.0 

957±16.2 

687±37.0 

5.42±0.59 



0%-10% 

327±9.5 

873±15.8 

624±32.4 

5.17±0.56 



10%-20% 

235±7.7 

597±13.4 

415±20.0 

4.28±0.47 



20%-30% 

166±6.3 

403±11.3 

274±15.1 

3.48±0.40 



30%-40% 

114T5.3 

263±10.1 

177±11.6 

2.74±0.34 



40%-50% 

75.0±4.5 

162±6.1 

110±9.2 

2.06±0.28 



50%-60% 

46.4±4.0 

91.5±6.2 

61.6±7.1 

1.38T0.23 



60%-70% 

26.1±3.5 

51.3±6.9 

31.6±5.0 

0.83T0.18 

Cu-bCu 

200 GeV 

0%-10% 

96.9±3.9 

238±12.2 

178±14.2 

3.00±0.36 



10%-20% 

74.3±3.9 

175±10.5 

123±9.9 

2.43±0.27 



20%-30% 

53.7±2.7 

121±8.7 

85.0±6.8 

2.00±0.25 



30%-40% 

39.9±3.8 

87.1±9.0 

57.7±4.6 

1.58±0.19 



40%-50% 

28.1±3.3 

59.0±7.9 

38.2±3.0 

1.24±0.17 

Au-I-Au 

62.4 GeV 

0%-10% 

317±6.1 

824±21.0 

405±32.4 

3.41±0.36 



10%-20% 

225±9.3 

560±17.4 

273±20.9 

2.95±0.30 



20%-40% 

131±8.5 

310±12.9 

151±13.1 

2.17±0.22 



40%-60% 

54.7±6.0 

118±8.0 

57.5±4.3 

1.31T0.13 

Cu-lCu 

62.4 GeV 

0%-10% 

95.9±2.1 

222±9.1 

122±8.9 

1.98±0.22 



10%-20% 

73.7±2.6 

164±8.4 

84.5±6.5 

1.65±0.19 



20%-30% 

55.2±2.5 

118±7.0 

58.0±4.5 

1.35±0.16 



30%-40% 

40.5±2.4 

83.6T6.7 

39.0±3.0 

1.10±0.13 



40%-50% 

28.2±2.2 

56.0±5.1 

25.5±2.0 

0.89±0.11 

Pb-fPb 

2.76 TeV 

0%-5% 

383±3.1 

1086±14.1 

1601±60 

11.5±1.43 



5%-10% 

330±4.6 

915±11.9 

1294±49 

10.5±1.27 



10%-20% 

261±4.4 

706±10.6 

966±37 

9.05±1.41 



20%-30% 

186±3.9 

488±8.3 

649±23 

7.35±1.21 



30%-40% 

129±3.3 

325±7.5 

426±15 

5.99±0.91 



40%-50% 

85.0±2.6 

205±5.9 

261±9 

4.69±0.75 



50%-60% 

52.8±2.0 

118±3.5 

149±6 

3.47±0.49 



60%-70% 

30.0T1.3 

60.9±2.0 

76±4 

2.11±0.35 



70%-80% 

15.8±0.6 

26.3±0.9 

35±2 

1.17±0.22 


A clear increase of the S'loss is seen in the 4-10 GeV/c 
region with the maximum being dependent on central¬ 
ity. Despite the ~10%-fold difference of y/s^ between 
RHIC and LHC, the trend is rather consistent, but more 
pronounced at the LHC and without a region of constant 
S'loss as is most evident in the PHENIX 0%-I0% data in 

Fig. 11 

The ALICE experiment recently published the spec¬ 
tra for charged pions for two centrality classes E3- We 
computed the fractional momentum loss for charged pi¬ 
ons and compared with those for charged hadrons as 
shown in Eig. H For peripheral collisions, we plot the 


results for charged hadrons in 60%-70% and 70%-80% 
bins. For 0%-5 % centrality, the Sioss for charged hadrons 
are systematically lower than that of charged pions at 
Pt <10 GeV/c, and both of them become similar above 
10 GeV/c. This observation is consistent with the en¬ 
hanced baryon production in px <10 GeV/c compared 
to mesons in the central collisions E3- Charged hadron 
spectra include protons, and thus the suppression is 
smaller for them in the medium px region. In the 60%- 
80% centrality, the charged pions and charged hadrons 
give similar results. This feature is again consistent with 
the observation of enhanced baryon production both at 
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FIG. 5. (Color online) dependence of Sioas for tt® in 
62 GeV Au+Au collisions using the spectra measured by 
PHENIX in 2010 Qj- SsysiTAA © pp norm) are Type-C errors 
and show the absolute amount that the data points would 
move. 



FIG. 7. (Color online) dependence of Sioss for charged 
hadrons in 2.76 TeV Pb+Pb collisions using the result from 
the ALICE experiment [3, [il. 4ys {Taa © pp norm) are 
Type-C errors and show the absolute amount that the data 
points would move. 




FIG. 6. (Color online) p^^ dependence of Sioss for 7r° in 
62.4 GeV Cu+Cu collisions using the spectra measured by 
PHENIX in 2005 [l3]. Ssys{TAA © pp norm) are Type-C er¬ 
rors and show the absolute amount that the data points would 
move. 


FIG. 8. (Color online) pf^ dependence of Sioss for charged 
pions in 2.76 TeV Pb+Pb collisions together with those for 
charged hadrons from the same collision system. The charged 
pion result is from the ALICE experiment 0. 


RHIC and LHC which only occurs in the central col¬ 
lisions. The ALICE experiment also published neutral 
pion data very recently, from which we calculated the 
5'ioss for the data set as shown in Figure El [Hi. 

The neutral pion results have finer centrality selec¬ 
tions, but have a limited pr range and larger uncertain¬ 
ties, therefore, they were not considered in further stud¬ 
ies of scaling variable dependence. We can see that the 
5'ioss for neutral pions are similar to that of charged pi¬ 
ons and hence are consistent with charged hadrons for 
Pt >10 GeV/c. 


B. Scaling variable dependence 

To understand how the fractional momentum loss 
changes with collision systems, we plot 5ioss against the 
scaling variables defined in the section |IT1 Figures [TOl 
and HD show the 5ioss as a function of Apart, fVqp, 
dNch/drj, and esj'To at = 7 and 12 GeV/c, respec¬ 
tively. Note that at these pf^ values, only data from 200 
GeV and 2.76 TeV are available. When a value at the ex¬ 
act was not available, we interpolated the fractional 
momentum loss from the closest two px points that we 
obtained in the previous section. The error bars repre- 
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FIG. 9. (Color online)py^ dependence of Sioss for neutral pions 
in 2.76 TeV Pb+Pb collisions using the result from the ALICE 
experiment [l^ . 


sent Type A and the boxes are Type B uncertainties; 
Type C uncertainties are not shown here. The scaling 
variable dependencies show clearer power-law behavior 
at pt = 12 GeV/c than at pp = 7 GeV/c, implying that 
the Sloss is dominated by a single source, i.e., hard scat¬ 
tering. At fixed the S\os& values for the Gu-I-Cu 

and Au-I-Au systems converge as IVpart grows. For the 
different values, a clear separation of S\oss values 

is seen even at the highest A^part, and the separation in¬ 
creases with increasing px (see Fig. [T^ . 

Figures [T2HIS1 show the same 5'ioss dependencies for ad¬ 
ditional pf^ values of 5-15 GeV/c. For the lowest two pi^ 
values, the results now also include Gu+Gu and Au+Au 
at = 62.4 GeV. Note that the PHENIX and AL- 

IGE data show parallel trends as a function of A^part, es¬ 
pecially at higher A^part- This fact, albeit the magnitudes 
are different, can be associated with the observation that 
ALICE and PHENIX data exhibit a similar A^part de¬ 
pendence of the dNch/dr]/{0.5NpaTt) shapes [l^. When 
looking at A"qp dependence, as expected from the discus¬ 
sion in the section explaining A^qp, the points are shifted 
up by a factor of 2-3 along the x-axis. The overall trends 
are similar as for A^part dependence, but the slopes are 
somewhat different. Comparing the data from different 
collision systems at the same ^/s^ reveals no significant 
improvement of the alignment from A^part to A^qp scaling. 
When we plot the 5'ioss against dNch/dp, the situation is 
different. 

At higher centralities (increasing dNch/dp) the LHC 
points line up very well with the 200 GeV RHIC Au+Au 
data, moreover, at higher pr the two results are con¬ 
sistent for all but the most peripheral collisions. This 
clearly shows that 5ioss scales with dNch/di], which is en¬ 
ergy density dependent and thus dependent. Fi¬ 

nally, plots of 5ioss as a function of eBjTblUlshow remark¬ 
able universal trends for the data from different systems 


from 200 GeV to 2.76 TeV. Among the scaling variables, 
dNch/dr] and esjA) seems to serve best across the colli¬ 
sion systems, especially between 200 GeV Au+Au and 
2.76 TeV collisions. This investigation shows that the 
5ioss does not scale with simple geometry descriptions 
across the , t>ut do scale with the quantities related 
to the energy density of the system, hence the opacity of 
the system is energy-density dependent. 

We have investigated 5ioss against the four scaling vari¬ 
ables at six p^ points including the two already shown in 
Figs. ITOl and ITT] The scaling plots at all are shown in 
Figs.[T2]-[ini For Pt of 5 and 6 GeV/c, we used the 2004 
data, because the 2007 data has a software threshold in 
Pt-, as mentioned earlier. At the same two lowest px, 
we also show the 5ioss scaling for 62.4 GeV Cu+Cu and 
Au+Au collisions. For higher px the 62.4 GeV points are 
not available owing to the lack of a p+p baseline. Devi¬ 
ations seen in the 62.4 GeV data may indicate that in 
the measured px range hard scattering is not completely 
dominant yet, in accordance with the observations of (7|. 

Lastly, to quantify the scaling trends, we fit S\oss for 
all four scaling variables and each collision system, except 
for ^/s^ = 62.4 GeV system, with a power-law function: 

Spx/PT = fdiSV/SV°r (7) 

where SV is one of the four scaling variables we used 
above, and the 5V° is the normalization factor intro¬ 
duced to cancel the dimension of the SV. We took the 
scaling variables for the most central LHC points as 5V°. 
Use of the power-law function is motivated by an energy 
loss model that predicts that AE/E oc Vpart^^^ [3l|. In 
the fitting process the statistical and systematic uncer¬ 
tainties were taken into account according to the pre¬ 
scription of [S^ . The errors on the scaling variable (hor¬ 
izontal errors in the plots) are not taken into account in 
the fitting, but they are small compared to the uncer¬ 
tainties of 5ioss values. 

The fit parameters a and /3 obtained by fitting 6px jpx 
vs Vpart and A(^p, plus dNch/dp and EBjA) to Eq. 0 
for Au+Au at = 200 GeV and Pb+Pb at 

= 2.76 TeV are shown in Fig. [161 All fit parameters, 
including for Cu+Cu, are tabulated in Table IVHl 

The fit parameters a and /3 are anti-correlated. At and 
above 10 GeV/c, the x^/nc?/ values become smaller and 
the powers a converge for all scaling variables, although 
they do not become fully consistent within uncertainties. 
Among the scaling variables, dNc^i/dr] is found to give 
relatively consistent a and /3 between two systems. The 
£BjTb, which is more related to the energy density of the 
system, also gives reasonably consistent numbers within 
uncertainties. More interestingly, £BjTo gives the a clos¬ 
est to 1.0 (linear scaling). The similarities are striking 
as is the fact that 5ioss obeys such a simple scaling with 
global observables over the entire px range where hard 
scattering is dominant. This implies that the empiri¬ 
cal fractional momentum loss and the assumed underly¬ 
ing energy loss of partons scale with energy density of 















12 


1 r 


10' r 


a. 

Q- H 

a 

CL 

CO 


(/) 


10-“' r, 


(a) S vs N ,,, 

' ' loss 

S,„ss@Pr = 7GeV/c 

, • *0 ® 
I *o° 


I ' 08 ° 


♦ Pb+Pb 2.76TeV, h*' 
PLB 720, 52 

OAu+Au 200GeV, It" 
PRL 101, 232301 

□Cu+Cu 200GeV, It" 
PRL 101,162301 


iJj td 


S|oss ''® 

S,os=@Pr = 7GeV/c 


I cPft 


t) 


(i)' 


10 


10 " 


10'’ 


10 ^ 10 


lO'' 


lO’’ 


part 


1 r 


10 ' r 


(b)S,„33VsN,, 
S,„ss@Pr = 7GeV/c 


, < 'o*® 
. I ' 0 ° 

>0 


(d) S, vs eR:T« 

' ’ loss 0 

S,o3s@Pr = 7GeV/c 


u- f 


,.o 


10 " 


dN^^dn 


PHENIX 


J Li. 




10 


10" 


10 '’ 


10 10 " 

eBi'to (GeV/fm^) 


FIG. 10. (Color online) Scaling variables dependence of Sioss at = 7 GeV/c. (a) shows Sioss vs A^part, (b) shows Sioss vs 
A^qp, (c) shows S'loss vs dNch/dp, and (d) shows S'loss vs £BjTo. A^qp are all calcnlated by PHENIX. Ssys{TAA © pp norm) is not 
shown in these plots. 
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FIG. 11. (Color online) Scaling variables dependence of Sioss at = 12 GeV/c. (a) shows Sioss vs Npart, (b) shows S'loss vs 
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the medium, independent of the collision energies or sys¬ 
tems, once is sufficiently high. We cross-checked 

our current result with one published earlier for a slightly 
different quantity [l^, and found consistent for 
= 200 GeV Au-I-Au collisions. 
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FIG. 12. (Color online) A^part dependence of the fractional momentum loss in bins of for various systems and 
Ssys{TAA © pp norm) is not shown in these plots. 



FIG. 13. (Color online) A^qp dependence of the fractional momentum loss in bins of pr for various systems and . Says{TAA 
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FIG. 14. (Color online) dNch/drj dependence of the fractional momentum loss. 5sys{TAA © pp norm) is not shown in these 
plots. 



FIG. 15. (Color online) eBjro dependence of the fractional momentum loss. 5ays{TAA © pp norm) is not shown in these plots. 
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IV. SUMMARY 

We have studied fractional momentum loss (iSioss 
=6pt/pt) over various systems and collision energies as 
a function of pt and four scaling variables: Apart, Aqp, 
dNch/dr] and esj'To- We found that the same universal 
function of dNch/dr] or Ebj^o describes ^loss at RHIC 
(\/^ = 200 GeV) and LHC = 2.76 TeV), while 

Apart and Aqp do not. This finding shows that the Sioss 
does not scale simply with system size across the 
but does scale with quantities related to the energy den¬ 
sity of the system, implying that the opacity of the sys¬ 
tem is energy-density dependent. We quantitatively eval¬ 
uated the slope of the universal curves for = 200 

and 2.76 TeV and again found that dNch/drj and eBj'Hj 
give relatively consistent a and (3 between two systems, 
and especially, that the the a for eBjT'o is close to 1.0 (lin¬ 
ear scaling). It is striking that S\oss obeys such a simple 
scaling with global observables over the entire pr range 
where hard scattering is dominant. This implies that the 
empirical fractional momentum loss and the assumed un¬ 
derlying energy loss of partons scale with energy density 
of the medium, independent of the collision energies or 
systems, once is sufficiently high. 

We propose that measurements of S'loss as well as the 
conventional Raa, in the future, would provide impor¬ 
tant additional information to investigate the global fea¬ 
ture of the energy loss of partons. 
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APPENDIX 

Tables of the centrality dependence of Spp /p^ and pa¬ 
rameters for fitting four different power-law functions for 
Au-|-Au and Cu-|-Cu data from the PHENIX experiment 
at RHIC and Pb-t-Pb data from the ALICE experiment 
at the LHC 
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TABLE IV. Centrality dependence of Spt/t^ in An+Au collisions at 200 GeV from 2007 and 2004 data from the 

PHENIX experiment at RHIC. 


Centrality 

P^T 

[GeV/c] 

2007 data 

5ptIp^ Stat error 

Syst error 

P^T 

[GeV/c] 

2004 data 

Spr/p^ Stat error 

Syst error 

0%-5% 


7.0 

0.216 

+ 0.004 

- 0.004 

+ 0.015 

- 0.013 


5.0 

0.202 

+ 0.003 

- 0.003 

+ 0.015 

- 0.013 



10.0 

0.209 

+ 0.005 

- 0.005 

+ 0.016 

- 0.014 


6.0 

0.206 

+ 0.004 

- 0.003 

+ 0.015 

- 0.013 



12.0 

0.204 

+ 0.007 

- 0.006 

+ 0.016 

- 0.013 


7.0 

0.216 

+ 0.002 

- 0.002 

+ 0.015 

- 0.013 



15.0 

0.157 

+ 0.012 

- 0.010 

+ 0.026 

- 0.021 






0%-10% 


7.0 

0.210 

+ 0.001 

- 0.001 

+ 0.016 

- 0.013 


5.0 

0.196 

+ 0.002 

- 0.002 

+ 0.015 

- 0.013 



10.0 

0.202 

+ 0.004 

- 0.004 

+ 0.016 

- 0.014 


6.0 

0.202 

+ 0.002 

- 0.002 

+ 0.015 

- 0.013 



12.0 

0.200 

+ 0.006 

- 0.005 

+ 0.016 

- 0.013 


7.0 

0.211 

+ 0.003 

- 0.003 

+ 0.015 

- 0.013 



15.0 

0.162 

+ 0.010 

- 0.009 

+ 0.026 

- 0.020 






10%-20% 


7.0 

0.172 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 


5.0 

0.165 

+ 0.002 

- 0.002 

+ 0.016 

- 0.014 



10.0 

0.162 

+ 0.005 

- 0.005 

+ 0.016 

- 0.014 


6.0 

0.171 

+ 0.002 

- 0.002 

+ 0.015 

- 0.013 



12.0 

0.168 

+ 0.007 

- 0.006 

+ 0.017 

- 0.014 


7.0 

0.180 

+ 0.003 

- 0.003 

+ 0.015 

- 0.013 



15.0 

0.128 

+ 0.012 

- 0.011 

+ 0.029 

- 0.022 






20%-30% 


7.0 

0.140 

+ 0.002 

- 0.002 

+ 0.017 

- 0.015 


5.0 

0.137 

+ 0.002 

- 0.002 

+ 0.016 

- 0.014 



10.0 

0.135 

+ 0.006 

- 0.005 

+ 0.016 

- 0.014 


6.0 

0.144 

+ 0.003 

- 0.002 

+ 0.016 

- 0.014 



12.0 

0.131 

+ 0.006 

- 0.006 

+ 0.019 

- 0.016 


7.0 

0.145 

+ 0.003 

- 0.003 

+ 0.016 

- 0.014 



15.0 

0.090 

+ 0.016 

- 0.014 

+ 0.034 

- 0.026 






30%-40% 


7.0 

0.110 

+ 0.002 

- 0.002 

+ 0.018 

- 0.015 


5.0 

0.120 

+ 0.002 

- 0.002 

+ 0.016 

- 0.014 



10.0 

0.108 

+ 0.006 

- 0.006 

+ 0.017 

- 0.015 


6.0 

0.122 

+ 0.003 

- 0.003 

+ 0.016 

- 0.014 



12.0 

0.113 

+ 0.007 

- 0.007 

+ 0.019 

- 0.016 


7.0 

0.126 

+ 0.004 

- 0.004 

+ 0.016 

- 0.014 



15.0 

0.071 

+ 0.020 

- 0.016 

+ 0.037 

- 0.027 






40%-50% 


7.0 

0.080 

+ 0.002 

- 0.002 

+ 0.018 

- 0.016 


5.0 

0.091 

+ 0.002 

- 0.002 

+ 0.017 

- 0.015 



10.0 

0.076 

+ 0.008 

- 0.007 

+ 0.017 

- 0.015 


6.0 

0.089 

+ 0.003 

- 0.003 

+ 0.017 

- 0.015 



12.0 

0.091 

+ 0.008 

- 0.007 

+ 0.020 

- 0.017 


7.0 

0.092 

+ 0.004 

- 0.004 

+ 0.017 

- 0.015 



15.0 

0.075 

+ 0.045 

- 0.027 

+ 0.037 

- 0.028 






50%-60% 


7.0 

0.055 

+ 0.003 

- 0.003 

+ 0.019 

- 0.016 


5.0 

0.062 

+ 0.003 

- 0.003 

+ 0.017 

- 0.015 



10.0 

0.056 

+ 0.010 

- 0.009 

+ 0.018 

- 0.016 


6.0 

0.064 

+ 0.004 

- 0.004 

+ 0.017 

- 0.015 



12.0 

0.064 

+ 0.011 

- 0.010 

+ 0.023 

- 0.019 


7.0 

0.072 

+ 0.005 

- 0.005 

+ 0.017 

- 0.015 



15.0 

0.029 

+ 0.027 

- 0.022 

+ 0.042 

- 0.031 






60%-70% 


7.0 

0.028 

+ 0.004 

- 0.004 

+ 0.019 

- 0.017 


5.0 

0.049 

+ 0.003 

- 0.003 

+ 0.017 

- 0.015 



10.0 

0.011 

+ 0.021 

- 0.019 

+ 0.028 

- 0.024 


6.0 

0.041 

+ 0.006 

- 0.005 

+ 0.018 

- 0.015 



12.0 

0.037 

+ 0.025 

- 0.022 

+ 0.046 

- 0.037 


7.0 

0.044 

+ 0.007 

- 0.006 

+ 0.018 

- 0.015 



15.0 

-0.098 

+ 0.046 

- 0.063 

+ 0.053 

- 0.077 
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TABLE V. Centrality dependence of Spr/p^ in Au+Au collisions at = 62.4 GeV and Cu+Cu collisions at at y/Sjyjy = 200 

and 62.4 GeV from the PHENIX experiment at RHIC. 


System 

Gentrality 

Pt^ 

[GeV/c] 


stat 

syst 

System 

Centrality 

PP 

Pt 

[GeV/c] 


stat 

syst 

\/^NN 



uncert. 

uncert. 




uncert. 

uncert. 

Au+Au 

0%-10% 

5.0 

0.115 

+ 0.010 

- 0.009 

+ 0.018 

- 0.015 

Cu+Cu 

30%-40% 

5.0 

0.034 

+ 0.002 

- 0.002 

+ 0.026 

- 0.021 

62.4 GeV 


6.0 

0.120 

+ 0.030 

- 0.023 

+ 0.019 

- 0.016 

200 GeV 


6.0 

0.033 

+ 0.004 

- 0.004 

+ 0.026 

- 0.021 







(continued) 


7.0 

0.036 

+ 0.007 

- 0.007 

+ 0.026 

- 0.021 


10%-20% 

5.0 

0.083 

+ 0.012 

- 0.010 

+ 0.019 

- 0.016 



10.0 

0.013 

+ 0.015 

- 0.013 

+ 0.031 

- 0.025 



6.0 

0.112 

+ 0.019 

- 0.016 

+ 0.019 

- 0.016 



12.0 

0.016 

+ 0.016 

- 0.015 

+ 0.035 

- 0.028 









15.0 

-0.001 

+ 0.028 

- 0.035 

+ 0.028 

- 0.035 


20%-40% 

5.0 

0.057 

+ 0.013 

- 0.012 

+ 0.020 

- 0.016 









6.0 

0.072 

+ 0.027 

- 0.021 

+ 0.020 

- 0.017 


40%-50% 

5.0 

0.015 

+ 0.004 

- 0.004 

+ 0.029 

- 0.024 









6.0 

0.022 

+ 0.006 

- 0.006 

+ 0.027 

- 0.022 

Cu+Cu 

0%-10% 

5.0 

0.102 

+ 0.001 

- 0.001 

+ 0.024 

- 0.020 



7.0 

-0.002 

+ 0.015 

- 0.016 

+ 0.034 

- 0.042 

200 GeV 


6.0 

0.103 

+ 0.002 

- 0.002 

+ 0.024 

- 0.020 



10.0 

0.033 

+ 0.021 

- 0.018 

+ 0.038 

- 0.031 



7.0 

0.098 

+ 0.004 

- 0.004 

+ 0.024 

- 0.020 



12.0 

0.023 

+ 0.029 

- 0.023 

+ 0.031 

- 0.025 



10.0 

0.074 

+ 0.008 

- 0.007 

+ 0.027 

- 0.022 



15.0 

0.060 

+ 0.099 

- 0.051 

+ 0.034 

- 0.026 



12.0 

0.076 

+ 0.009 

- 0.008 

+ 0.027 

- 0.022 









15.0 

0.062 

+ 0.020 

- 0.017 

+ 0.029 

- 0.023 

Cu+Cu 

0%-10% 

5.0 

0.041 

+ 0.012 

- 0.010 

+ 0.028 

- 0.022 


10%-20% 

5.0 

0.078 

+ 0.002 

- 0.002 

+ 0.024 

- 0.020 

62.4 GeV 


6.0 

0.057 

+ 0.034 

- 0.025 

+ 0.030 

- 0.023 



6.0 

0.077 

+ 0.003 

- 0.003 

+ 0.024 

- 0.020 









7.0 

0.075 

+ 0.005 

- 0.004 

+ 0.025 

- 0.020 


10%-20% 

5.0 

0.036 

+ 0.013 

- 0.011 

+ 0.027 

- 0.021 



10.0 

0.054 

+ 0.009 

- 0.008 

+ 0.028 

- 0.022 



6.0 

0.048 

+ 0.035 

- 0.026 

+ 0.030 

- 0.023 



12.0 

0.065 

+ 0.010 

- 0.009 

+ 0.028 

- 0.022 









15.0 

0.011 

+ 0.025 

- 0.021 

+ 0.036 

- 0.027 


20%-30% 

5.0 

0.016 

+ 0.015 

- 0.013 

+ 0.029 

- 0.022 









6.0 

0.024 

+ 0.031 

- 0.024 

+ 0.028 

- 0.022 


20%-30% 

5.0 

0.051 

+ 0.002 

- 0.002 

+ 0.025 

- 0.021 









6.0 

0.054 

+ 0.004 

- 0.004 

+ 0.025 

- 0.021 


30%-40% 

5.0 

0.005 

+ 0.028 

- 0.024 

+ 0.056 

- 0.044 



7.0 

0.048 

+ 0.006 

- 0.006 

+ 0.026 

- 0.021 



6.0 

-0.010 

+ 0.127 

- 0.163 

+ 0.137 

- 0.180 



10.0 

0.028 

+ 0.011 

- 0.010 

+ 0.029 

- 0.023 









12.0 

0.055 

+ 0.014 

- 0.012 

+ 0.029 

- 0.023 


40%-50% 

5.0 

-0.019 

+ 0.018 

- 0.021 

+ 0.026 

- 0.033 



15.0 

0.034 

+ 0.028 

- 0.022 

+ 0.029 

- 0.023 



6.0 

-0.034 

+ 0.035 

- 0.050 

+ 0.019 

- 0.024 
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TABLE VI. Centrality dependence of 6pt/p^ Pb+Pb collisions at = 2.76 TeV from the spectra measured by the ALICE 

experiment at the LHC [l^. [itI. [3^ . 


Centrality 


[GeV/c] 

SpTlp^r^ 

Stat error 

Syst error 

Centrality 


[GeV/c] 

+ t / p ?? 

Stat error 

Syst error 

0%-5% 


5.0 

0.241 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 

30%-40% 


10.0 

0.187 

+ 0.002 

- 0.002 

+ 0.018 

- 0.016 



6.0 

0.270 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 

(continued) 


12.0 

0.173 

+ 0.002 

- 0.002 

+ 0.018 

- 0.016 



7.0 

0.293 

+ 0.001 

- 0.001 

+ 0.015 

- 0.014 



15.0 

0.154 

+ 0.004 

- 0.004 

+ 0.020 

- 0.017 



10.0 

0.316 

+ 0.001 

- 0.001 

+ 0.015 

- 0.013 









12.0 

0.303 

+ 0.001 

- 0.001 

+ 0.015 

- 0.013 

40%-50% 


5.0 

0.141 

+ 0.001 

- 0.001 

+ 0.019 

- 0.017 



15.0 

0.282 

+ 0.002 

- 0.002 

+ 0.016 

- 0.014 



6.0 

0.153 

+ 0.001 

- 0.001 

+ 0.019 

- 0.016 









7.0 

0.158 

+ 0.001 

- 0.001 

+ 0.019 

- 0.016 

5%-10% 


5.0 

0.229 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 



10.0 

0.148 

+ 0.002 

- 0.002 

+ 0.019 

- 0.017 



6.0 

0.255 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 



12.0 

0.142 

+ 0.003 

- 0.003 

+ 0.019 

- 0.017 



7.0 

0.277 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 



15.0 

0.123 

+ 0.005 

- 0.005 

+ 0.021 

- 0.018 



10.0 

0.293 

+ 0.001 

- 0.001 

+ 0.015 

- 0.014 









12.0 

0.281 

+ 0.002 

- 0.002 

+ 0.016 

- 0.014 

50%-60% 


5.0 

0.116 

+ 0.001 

- 0.001 

+ 0.020 

- 0.017 



15.0 

0.259 

+ 0.003 

- 0.002 

+ 0.017 

- 0.015 



6.0 

0.122 

+ 0.001 

- 0.001 

+ 0.020 

- 0.017 









7.0 

0.130 

+ 0.002 

- 0.002 

+ 0.020 

- 0.017 

10%-20% 


5.0 

0.211 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 



10.0 

0.118 

+ 0.003 

- 0.003 

+ 0.020 

- 0.018 



6.0 

0.236 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 



12.0 

0.105 

+ 0.004 

- 0.004 

+ 0.021 

- 0.018 



7.0 

0.253 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 



15.0 

0.084 

+ 0.007 

- 0.007 

+ 0.022 

- 0.019 



10.0 

0.263 

+ 0.001 

- 0.001 

+ 0.016 

- 0.014 









12.0 

0.252 

+ 0.002 

- 0.001 

+ 0.016 

- 0.014 

60%-70% 


5.0 

0.091 

+ 0.002 

- 0.002 

+ 0.021 

- 0.019 



15.0 

0.228 

+ 0.002 

- 0.002 

+ 0.018 

- 0.015 



6.0 

0.094 

+ 0.002 

- 0.002 

+ 0.021 

- 0.019 









7.0 

0.094 

+ 0.003 

- 0.002 

+ 0.021 

- 0.019 

20%-30% 


5.0 

0.190 

+ 0.001 

- 0.001 

+ 0.018 

- 0.015 



10.0 

0.086 

+ 0.004 

- 0.004 

+ 0.022 

- 0.019 



6.0 

0.210 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 



12.0 

0.080 

+ 0.006 

- 0.006 

+ 0.022 

- 0.019 



7.0 

0.224 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 



15.0 

0.071 

+ 0.011 

- 0.010 

+ 0.023 

- 0.020 



10.0 

0.224 

+ 0.001 

- 0.001 

+ 0.017 

- 0.015 









12.0 

0.212 

+ 0.002 

- 0.002 

+ 0.017 

- 0.015 

70%-80% 


5.0 

0.075 

+ 0.003 

- 0.003 

+ 0.023 

- 0.020 



15.0 

0.190 

+ 0.003 

- 0.003 

+ 0.019 

- 0.016 



6.0 

0.074 

+ 0.003 

- 0.003 

+ 0.024 

- 0.020 









7.0 

0.077 

+ 0.004 

- 0.004 

+ 0.023 

- 0.020 

30%-40% 


5.0 

0.168 

+ 0.001 

- 0.001 

+ 0.018 

- 0.016 



10.0 

0.068 

+ 0.006 

- 0.006 

+ 0.024 

- 0.021 



6.0 

0.183 

+ 0.001 

- 0.001 

+ 0.018 

- 0.016 



12.0 

0.081 

+ 0.010 

- 0.009 

+ 0.024 

- 0.020 



7.0 

0.195 

+ 0.001 

- 0.001 

+ 0.018 

- 0.015 



15.0 

0.054 

+ 0.019 

- 0.017 

+ 0.026 

- 0.022 
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TABLE VII. Parameters from fitting the indicated power-law functions for Spr/pr to the data as a function of p^ for Au+Au 
collisions from 2004 and 2007 data and for Cu-fCu collisions from 2005 data at = 200 GeV. 


System 

■\/^NN 

year 

hadron 

Spt/pt = 

P^T 

a 

P 



X^/ndf 

Au-I-Au 

200 GeV 

2004 


/ 3 (Arpart/Ar°,rt)“ 

5 GeV/c 

0.529l°°ll 

2 14+004 
J-+_0.03 

X 

10-1 

25.45/5 






6 GeV/c 

0.543l°°l® 

9 90+0.04 
Z.ZO_o_o 4 

X 

10-1 

15.56/5 






7 GeV/c 

U. 0 ^O_Q 020 

0 09+0.05 
Z.OZ_o_o 4 

X 

10-1 

7.11/5 






5 GeV/c 

U.^OO_o 010 

9 1 O+0.04 

4 .io_o.o 4 

X 

10-1 

23.35/5 






6 GeV/c 

0.475l°°« 

0 Q-y+0.04 
-0.04 

X 

10-1 

15.23/5 






7 GeV/c 

0.480l°°l^ 

9 qf:’+0.05 

z.ou_o_o 5 

X 

10-1 

7.50/5 





PidN,k/dp/dN°Jdpr 

5 GeV/c 

0.4451°°°® 

Q 01 +0 06 

O.Ui_o.o6 

X 

10-1 

27.78/5 






6 GeV/c 

0.456l°°« 

q 1 r+0.08 
o.-Ltj_o 07 

X 

10-1 

18.56/5 






7 GeV/c 

0.460l°°l^ 

Q QQ+O IO 
O.OU_o.o 9 

X 

10-1 

8.50/5 





/ 3 (ero/e°ro)“ 

5 GeV/c 

U.OiO_o.oi8 

q yq+0.09 

0.1 l>_o.09 

X 

10-1 

14.67/5 






6 GeV/c 

U.OOZ_0.025 

4.00l°;l® 

X 

10-1 

3.79/5 






7 GeV/c 

U.OOlt_0 032 

4.17l°;l° 

X 

10-1 

4.23/5 

Au-I-Au 

200 GeV 

2007 

7r° 


10 GeV/c 

U.UOZ_0.035 

0 OQ+O.Ut) 

Z.ZO_o.o6 

X 

10-1 

3.31/5 






12 GeV/c 

0.56ll°°« 

9 10+0.07 

^•-Ly_o.o 7 

X 

10-1 

1.75/5 






15 GeV/c 

0.795l°:lfi 

i.oO_0.i3 

X 

10-1 

4.68/5 





m^p/Kr 

10 GeV/c 

A c i::9+0.032 
U.OOZ_o 031 

9 90+O.O6 

z.zo_0_06 

X 

10-1 

3.32/5 






12 GeV/c 

0.490l°°°® 

0 00+0.07 

z.zz_o_o7 

X 

10-1 

1.78/5 






15 GeV/c 

0.695l°;lg 

i.yu_o 14 

X 

10-1 

4.74/5 





PidN,k/dp/dN°Jdpr 

10 GeV/c 

0 cioq+0-030 

U.OZO_o.o29 

q qq + 0.15 

O.OO_0.i4 

X 

10-1 

3.72/5 






12 GeV/c 

0 471+0 033 
'J.lt/ -l-_o.032 

q 1 q+0.17 
O.i0_o 15 

X 

10-1 

1.59/5 






15 GeV/c 

A +0-124 

0.661Io,ii7 

O.UO_o.42 

X 

10-1 

4.69/5 





/ 3 (ero/e°ro)“ 

10 GeV/c 

1.0201°°°° 

4.541°-®° 

X 

10-1 

2.05/5 






12 GeV/c 

0.8921°°°^ 

4.051°:®° 

X 

10-1 

2.43/5 






15 GeV/c 

1.3001°;®°® 

^•O 0 _o. 9 i 

X 

10-1 

4.36/5 

Cu-fCu 

200 GeV 

2005 

7r° 

l3{Np^rt/K^,,r 

5 GeV/c 

1 01 A+O.U 4 t) 

i.ZlU_o.045 

c 4 c:+ 0-41 
^•4:0_o 37 

X 

10-1 

8.28/3 






6 GeV/c 

1 1 or>+0.082 

i.ioU_0 079 

r 01+0.70 

^•^-L_o.60 

X 

10-1 

1.48/3 






7 GeV/c 

i.ZUU_o.i 4 i 

r 1 r^+1.31 
0--*-/_i 01 

X 

10-1 

2.92/3 





m^p/Kr 

5 GeV/c 

i.UOU_o 039 

tx. 91 +0.38 

O.Zl_o 35 

X 

10-1 

9.71/3 






6 GeV/c 

i.UOU_o 069 

4 QQ+0-65 

^•yy_o.56 

X 

10-1 

1.69/3 






7 GeV/c 

i.UOU_o 124 

4 Q4+1-21 

^•y4:_o 94 

X 

10-1 

3.07/3 





PidN,k/dp/dN°Jdpr 

5 GeV/c 

U.y4:U_o 035 

0 99 + 0.74 

O.zz_o 67 

X 

10-1 

15.46/3 






6 GeV/c 

0 Q1 7+0 063 
u.yj- 1 _0.061 

7.80ll®° 

X 

10-1 

2.26/3 






7 GeV/c 

0 Q31+0113 

u.yoj._o 108 

7 70+2.39 

/./U_i 77 

X 

10-1 

3.81/3 





/ 3 (ero/e°ro)“ 

5 GeV/c 

1 (^7Q+0-063 

i.O/U_o 061 

9.831°-°° 

X 

10-1 

15.29/3 






6 GeV/c 

1.630l°:“® 

9.28ll°t 

X 

10-1 

1.92/3 






7 GeV/c 

1.6501°-®°® 

9.18lli°5 

X 

10-1 

3.88/3 
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TABLE VIII. Parameters from fitting the indicated power-law functions for 5pt/pt to the data as a function of p^ for Pb-|-Pb 
collisions at =2.76 TeV. 


System 

year 

hadron 

SptIpt = 

P^T 

a 

p 



X^/ndf 

Pb-fPb 2.76 TeV 

2010-11 


/3(Afpart/fVpVt)“ 

5 GeV/c 

n oc;7+0-004 

u.o<Ji _o.o04 

^•^^-0.04 

X 

10-1 

44.19/7 





6 GeV/c 

n 070+0.004 

u.o/o_Q QQ3 

2.741^4 

X 

10-1 

90.44/7 





7 GeV/c 

0.3981Hm 

2.961^4 

X 

10-1 

70.86/7 





10 GeV/c 

0.490tH2^ 

'J.-LU_0.04 

X 

10-1 

10.32/7 





12 GeV/c 

0 c;r)7+0 008 
U.JU/ _ 0.008 

o.U4:_0_04 

X 

10-1 

11.41/7 





15 GeV/c 

KJ.OOl _o.014 

Q 09 + 0.05 
Z.OZ_0.04 

X 

10-1 

2.29/7 




m^p/K^r 

5 GeV/c 

U.OZU_o.o03 

9 44+^ '^^ 
^•^^-0.04 

X 

10-1 

34.51/7 





6 GeV/c 

0 ‘^‘^Q+0 003 

U.OOy_Q QQ3 

9 7Q+0.04 
' 'J-0.04 

X 

10-1 

71.06/7 





7 GeV/c 

0 Q(^Q+0 003 

u.oao_o.oo3 

9 Qc;+0.05 
z.yj_0.04 

X 

10-1 

59.14/7 





10 GeV/c 

U.4:^U_q 005 

'J.J-U_0.04 

X 

10-1 

9.62/7 





12 GeV/c 

U . 4 : OO _ 0.007 

o.U+_0.04 

X 

10-1 

13.94/7 





15 GeV/c 

O.SOllHls 

9 qo+0.05 
Z.OO_0.04 

X 

10-1 

2.30/7 




/3{dN,H/dp/dN°Jdpr 

5 GeV/c 

0.298tg;°S3 

^•+O_0.04 

X 

10-1 

66.71/7 





6 GeV/c 

n *^10+0.003 
U.OJ.O_0.003 

9 77+0.04 

' ' -0.04 

X 

10-1 

145.00/7 





7 GeV/c 

0.329tHo3 

2.981^5 

X 

10-1 

123.28/7 





10 GeV/c 

0 404+0-0*^^ 
U.4:U4:_o,005 

0 iq+ 0.05 
'j.j-y_o.o4 

X 

10-1 

30.94/7 





12 GeV/c 

0.417+nS6 

O.UU_0.04 

X 

10-1 

26.21/7 





15 GeV/c 

0.455tHn 

9 Qc :+ 0-05 
Z.OO_o 04 

X 

10-1 

5.76/7 




/3(eTo/e°ro)“ 

5 GeV/c 

0 f ^ 7 f ;+0 006 

U.O/O_0.006 

2 ^ o + 0.04 
^•^ 9 _ o . o 4 

X 

10-1 

53.83/7 





6 GeV/c 

U.O14:_0,006 

9 7Q+0.04 
' ’J-0.04 

X 

10-1 

91.36/7 





7 GeV/c 

U.O4:y_0.006 

^.yu_0.04 

X 

10-1 

79.47/7 





10 GeV/c 

0 7 QQ +0 009 
u. I ^^_ 0.009 

q 1 7+0.05 

O.J-1 _o.04 

X 

10-1 

32.58/7 





12 GeV/c 

u.oz,y_0.0i3 

O.UO_0.04 

X 

10-1 

30.78/7 





15 GeV/c 

u.yuy_0.023 

9 Qq + 0.05 
Z.OO_o.o4 

X 

10-1 

6.28/7 
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